Quantum Electrodynamical Photon Splitting in Magnetized Nonlinear Pair Plasmas 
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We present for the first time the nonlinear dynamics of quantum electrodynamic (QED) pho- 
ton splitting in a strongly magnetized electron-positron (pair) plasma. By using a QED corrected 
Maxwell equation, we derive a set of equations that exhibit nonlinear couplings between electro- 
magnetic (EM) waves due to nonlinear plasma currents and QED polarization and magnetization 
effects. Numerical analyses of our coupled nonlinear EM wave equations reveal the possibility of 
a more efficient decay channel, as well as new features of energy exchange among the three EM 
modes that are nonlinearly interacting in magnetized pair plasmas. Possible applications of our 
investigation to astrophysical settings, such as magnetars, are pointed out. 

PACS numbers: 52.35.Mw, 94.20.wf 



Recently, there has been a great deal of interest [ij in 
the investigation of effects associated with radiation pres- 
sure and quantum vacuum fluctuations in nonhnear me- 
dia. Such studies are of importance in astrophysical en- 
vironments, where copious amounts of electron-positron 
pairs exist due to numerous physical processes Elas- 
tic photon-photon scattering is traditionally described 
within quantum electrodynamic (QED) [l|, However, 
observable effects of elastic photon-photon scattering 
among real photons have so far not been reported in the 
laboratory [l|, Q . For astrophysical systems 0, 0] the sit- 
uation is different, since the large magnetic field strength 
in pulsar and magnetar Q environments changes the dia- 
magnetic properties of vacuum significantly [7[ , and leads 
to phenomena such as frequency down-shifting The 
latter is a result of photon splitting [1, Q , and the pro- 
cess may be responsible for the radio silence of magnetars 
0, Moreover, the propagation of electromagnetic 

waves in a relativistically dense electron gas [llj and in 
a relativistic electron-positron gas have been dis- 
cussed, leading to the important effect of gamma photon 
capture and pair plasma suppression around pulsars [ist . 

In this Letter, we present the nonlinear photon split- 
ting of electromagnetic (EM) waves propagating perpen- 
dicularly to a strong external magnetic field Bq in a pair 
plasma. Due to the QED effect [iJ, a photon in vacimm 
can decay into a backscattered and a forward scattered 
photon, where the latter two photons have polarizations 
perpendicular to that of the original photon [1,0] . Noting 
that significant pair-production (l^ occurs in astrophys- 
ical settings (viz. in pulsar and magnetar environments), 
we demonstrate here a novel possibility of a nonlinear 
decay interaction, due to a competition between QED 
and plasma nonlinearities. We note that most of pre- 
vious investigations 15, [lB|, including both QED and 



plasma effects, have been limited to linear EM wave prop- 
agation. Here we derive three dynamical equations with 
nonlinear couplings between photons with different polar- 
izations. From these coupled mode equations, the QED 



cross-section for photon-splitting [8| can be deduced in 
the limit of zero plasma density. We discuss applications 
of our results to magnetar atmospheres. 

Photon-photon scattering can be described by the 
Heisenberg-Euler Lagrangian [13, L = eo{E'^ — 



2^2)2 + 7(cE • B)2]. Here 



(a/907r)(l/eo£'2j.jj, a = e'^/^TTCohc is the fine-structure 
constant, i?crit = m^c^ /eh ~ 10^^ V/m is the critical field 
[l[, /i is the Planck constant, m is the electron mass, eo 
is the vacuum permittivity, and c is the vacuum speed 
of light. The last two terms in the Heisenberg-Euler 
Lagrangian represent the effects of the vacuum polariza- 
tion and magnetization. The QED corrected Maxwell 
equations can then be written in their classical form us- 
ing D = eoE + P and H = c^eoB - M, where [H 



- 9f2 



2b2)E + 7c2(E • B)B] and M 



2c^elK[-2{E^ - c^B^)B + 7(E • B)E], which are valid for 
|E|, c|B| < i:;crit and u ^ Ue = mc^/h « 8 x 10^° rad/s. 

Next, we study wave propagation perpendicular to an 
external magnetic field Bq = BqZ in an electron-positron 
plasma, letting all variables depend on {x,t). Assuming 
that the charge density is negligible, the wave equation 
for the electric field E then reads 

{df - c'dl) E = [da c^WlP.] , (1) 

where j 9tP + V x M + Y.e,p 9"^, |E|/c, |B| < Bq, v 
denotes the average (fluid) velocity, and the sum is over 
the electron and positron contributions. The latter are 
determined from the relativistic equation of motion 



(ft -f V V) (7v) = (g/m) (E + v x B) 



(2) 



We assume that the EM wave frequency uj and the elec- 
tron (positron) plasma frequency Wpe(p) are much smaller 
than the magnitude of the electron (or positron) gy- 
rofrequency |ti-'ce(p)| = = eBo/m, relevant for pulsar 
and magnetar atmospheres (2^ . This ordering will make 
charge density oscillations negligible, as the longitudinal 
motion for both the electrons and positrons will be given 



2 



by the E X B-drift to leading order. Next, we linearize 
and Fourier decompose Eqs. ^ and ([2]) to obtain the 
dispersion relations for the EM waves propagating per- 
pendicular to z [i^l 

Lu^ « fc^c^ (1 - 80 , (3a) 
Lu'' ^k^c^il-UO+^l (3b) 

where ujp ~ {ujp^ + w^^)^/^ is the plasma frequency of the 
pair plasma. Here, ^ = neoc^BQ ^ 1, and we assume 
that ujp ^ We have omitted the contribution 

proportional to u'^Up/tu'^ in (j3ap . which is smaller than 
the plasma contribution proportional to uj^ in (j3bp . We 
note that the EM mode described by ((3a|) has the elec- 
tric field perpendicular to Bq in the y-direction (approxi- 
mately) , whereas the EM mode described by ()3bp has the 
electric field parallel to Bq . Henceforth, the two different 
polarizations will be denoted by the subscripts _L and ||, 
respectively. In [l^ and the linear effects from the 
combined QED and plasma effects are discussed in detail. 

Next, we represent the EM waves as E{x, t) exp{ikx — 
iujt) + complex conjugate, and the slowly varying am- 
plitudes are denoted by tilde. Our aim is to investi- 
gate photon splitting [1, 0] , a parametric process where 
one photon with perpendicular polarization decays into 
two photons with parallel polarizations. Denoting the 
latter waves with indices 1 and 2, the energy and mo- 
mentum conservation relations (matching conditions) are 
w_L = +^^^211 and fc_L = fci|| -I- fc2||. We point out that 
in addition to this process, QED allows a decay of the 
type = + ^2\\ 0- However, the scattering am- 
plitude of this process is suppressed by a factor of the 
order (2] ■ The presence of a plasma may in principle 
change this ordering and also add new decay channels, 
but in the strongly magnetized high-frequency regime 
considered here, uj,lOc S> ujp, we note that this is not 
the case. For Wp = 0, the simultaneous fulfilment of the 
matching conditions and the dispersions relations ((3a|) 
and (|3b[) requires that one of the EM waves is backscat- 
tered, i.e. either < or k2\\ < 0. Furthermore, 
^ -C 1 means that the backscattered EM wave has a 
much smaller frequency than lu_i_. On the other hand, for 
Up ^ 0, the matching conditions also allow for both the 
decay products to be scattered in the forward direction. 
Next, we divide all quantities into unperturbed and per- 
turbed quantities, i.e. B = Bo + Bi where the perturbed 
part fulfills |Bi| ^ |Bo|, and similarly for the polariza- 
tion, magnetization and density. We then include the 
resonant second order nonlinear terms from P and M in 
Maxwell's Eqs. (noting that P = Po-l-Pi-f P2, where P2 
is second order in the perturbed EM-field, etc.) together 
with the second order terms in ([2]) (i.e. the Lorentz force 
and the convective derivative) , and the nonlinear part of 
the current density in the right-hand side of ([T]). From 
the continuity equation dtpi = —podxVx we solve for the 
density. After straightforward algebra, where the linear 



expressions are substituted into the nonlinear terms, we 
obtain our coupled mode equations 

dtE^ + Vg^dxE^ = u^CE^E2\\/E,,,,, (4a) 
dtE^ + = u^CE;^^E^/E„,,, (4b) 

dtE^w + Vg2dxE2\\ = uj2\\CE*^^E^/E„iu (4c) 

where Vgj = dujj/dkj is the group speed (j equals 
_L, 1|| and 2||), and the asterisk denotes complex con- 
jugate. The coupling strength is C = Cpi -I- Cqed, 

where Cpi = [{a/dOnO^^^ {kj_c/ujj){ujp/uJi\\uj2\\) is 
due to the plasma nonlinearities [25[ and Cqed = 

2i (ae/90^)'/' [I0(,k^c/iu^) + 7 (fci||cM|| + fc2||cM||)] 
is due to QED nonlinear interactions [26|. During 
certain conditions and for sufficiently long times, the 
multi-scale expansion behind Eqs. Q can break down 
due to the growth of higher order terms. However, we 
will assume that, due to convective stabilization (see 
below), such effects are negligible and the evolution is 
well described by ([4]). 

Let us now discuss the relative importance of the QED 
and plasma effects in various regimes. 

Case 1: In the regime ujp/uj±^ <^ ^, the matching 
conditions, the dispersion relations, and ^ <C 1 give 
^2|| ~ —3^k±/2 and uj2\\ ~ 3^a;^/2, where the last ap- 
proximation is valid to first order in ^, and we have cho- 
sen uj2\\ as the low- frequency backscattered wave with 
k2\\ < 0. We note that the condition ujp/uj± ^ £, ensures 
that the QED effect dominates over the plasma effects. 

Case 2: Increasing the plasma density starts to affect 
the linear properties of the low- frequency EM mode first. 
However, we note that the expression UJ2\\ = i£,u!±/2 even 
holds when the plasma effect dominates over the QED ef- 
fect in the dispersion relations for the low-frequency EM 
mode. Increasing the plasma frequency to the regime 
ujp/uj± ^ ^, the nonlinear coefficients Cqed and Cpi be- 
come comparable and the matching conditions further 
imply that 0^211 ~ ~ £^uj± [see the note after Eqs. (3)]. 
Furthermore, we note that for cOp > 3^w^ we would have 
forward scattering instead of backscattering of the low- 
frequency EM mode. 

Case 3: For larger plasma frequency {ujp/uj± ^ ^) the 
QED-contribution to the frequency uj± may still domi- 
nate over the plasma contribution, but |Cqed| ^ |Cpi|, 
and for the nonlinear wave interaction we can thus omit 
the QED effect in this regime. 

We note that the nonlinear system ^ has the con- 
served energy integral £ ^ J {\E±\'^ + {E^ p -I- \E2\\ P) dx 
when = uji\\ + uj2\\, and that the two other linearly 
independent constants of motion are A/i = / {\E±\'^/uj± — 
\E,\\\yLu,\\)dx and Ar2 = Ji\E^\Vu;^ - |£;2|| l'M||) ^a;, 
corresponding to the Manley-Rowe relations. The con- 
stants of motion are used as a check of the numeri- 
cal calculations presented below. We first make a lin- 
ear stability analysis in the presence of a pump EM 



3 




7? 



Mi 



-2.5 
2 



10.05 1 
0.045 0.9 
0.04 0.8 
0.035 0.7 
■0.03 0.6 
0.025 0.5 
0.02 0.4 
0.015 0.3 
0.01 0.2 
0.005 0.1 



1 

10,045 0.9 

0.04 0.8 

0.035 0.7 

0.03 0.6 

0.025 0.5 

0.02 0.4 

0.015 0.3 

0.01 0.2 

0.005 0.1 



0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 



-0.5 0.5 
xk,/10^ 



-1 -0.5 0.5 
xk,/10^ 



-1 -0.5 0.5 
xk/IO^ 



FIG. 1: The decay of the pump mode E± — E±/EctH 
into a forward scattered mode = E±/EcTit and a back- 
scattered mode -E211 = Ej_/Ecrit- InitiaUy, the pump is set 
to E± — 0.02, while Ei and E2 is set to a low-level random 
noise. After the initial exponential decay, the energy is trans- 
ferred between the pump and the two EM sidebands in a semi- 
periodic and chaotic manner. We used uip = and ^ = 0.01, 
yielding fci|| = 1.017fc±, fcan = -0.017fc_L, uj± = 0.959fcxc, 
cjill — 0.943fcxc, and uj2\\ = O.OlOfcxc. 
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FIG. 3: The decay of the pump mode E± = E±/EcTit in a 
strongly magnetized plasma. The pump was initially set to 
a localized pulse E± = 0.05 exp[-(a:A:x - 10^)V2.5 x 10^]. 
The pump decays into a forward scattered sideband Siy 
and a sideband E2\\, which has almost zero group speed. 
We used ^ = 0.01 and ujp — 3£,io± — 0.03a;x, yielding 
fci|| = k±, k2\\ ^ 0, u!± = 0.959fcxc, cjiy = 0.929fcxc, and 
L02\\ — iUp ^ 0.030A;xC. 



main —10^ < xk± < 10^ with periodic boundary condi- 
tions, and 20 000 steps to advance the solution in time. A 
pseudo-spectral method is used to approximate the spa- 
tial derivatives and a 4*'^-order Runge-Kutta method for 
the time stepping. In Figs. 1 and 2, the plasma is ab- 
sent, so that the only EM wave couplings are due to the 
QED effect. We used C = 0.01, yielding fci|| = 1.017A:j., 
k2\\ = -0.017fc_L, ujj_ = 0.959fc_LC, wiy = 0M3k±c, 
and CJ211 = G.016A:^c. In Fig. 1, we present the evo- 
lution of an initially homogeneous beam of amplitude 
E± = 0.02i?crit- Initially, the two daughter waves grow 
exponentially, followed by a nonlinear oscillatory phase. 
Figure 2 exhibits the nonlinear dynamics of a localized 
wave packet. We observe the decay of the EM pulse into 
a forward scattered wave Eiu and a backscattered wave 



FIG. 2: The decay of the pump mode E± — E±/EcTit- 
The pump was initially set to a localized pulse E± — 
0.05exp[-(a:ifcx - 10^)V2.5 x lO'']. The pump decays into a E2\\. In Fig. 3, we show the dynamics of a localized EM 
forward scattered sideband Siy and aback-scattered sideband 
£211- We used ujp = and ^ = 0.01, yielding fciy = 1.017fcx, 
/tail = -0.017itx, cj± = 0.959fcxc, = 0.943fcxc, and 

LO2H = 0.016A;xc. 



wave. Thus, we consider the decay of a homogeneous 
intense wave Ej_ = E±o, where \E±q\ > I-B2I1I, 
into daughter EM waves = cxp{iKx — i^lt) 

and -E'211 = E2\\exp{—iKx + int). From (g]), we ob- 
tain the nonlinear dispersion relation (fi — Vg2\\K){Q — 
Vgi\iK) = -uJ^UJ2\\\C\'^\EJ_o\'^/E^^■^. From the latter, the 
growth rate of the daughter EM waves is obtained as 
r = [u.i||W2|l|Cni?^olV^^cnt-Ki|| +«s2|l)'^V4]^/' for 
wavenumbers iC < 2(wi||a;2||)^^^|C-B_Lo|/|'!;g2||+^^gi||l-E^crit- 
In order to study the dynamics of an intense electro- 
magnetic beam in a strong magnetic field, we numerically 
solve the coupled Eqs. and display the results in Figs. 
1-3. Wc use 1000 grid points to resolve the numerical do- 



pulse when the plasma effect is important. We consider 
the particular case ujp = 3^a;_L, so that the low- frequency 
daughter EM wave E2\\ has approximately zero group 
speed and a frequency that equals the plasma frequency. 
We thus used = 0.01 and ujp — 3^uj± ~ 0.03Li;_L, yield- 
ing A:i|[ = k±, k2\\ ~ 0, uj± = 0.959A:^c, = 0.929A;j_c, 
and LU2\\ = i-Up 0.030fcj^c. The energy of the pump 
E±_ is transferred to a forward scattered wave and 
zero-group speed waves E2\\- 

The present investigation is of relevance for EM wave 
propagation in the vicinity of pulsars and magnetars. For 
example, the radio silence of magnetars is assumed to be 
connected with the photon-splitting in the strong magne- 
tar fields (10^ - 10" T) [i,[ia|. Photon splitting can sup- 
press the creation of electron-positron pairs [10| , but we 
still expect the presence of an electron-positron plasma 
14| in such environments. The Goldreich- Julian den- 
sity is given hy ^ naj = 7 x 10^^ {0. 1 /t) {Bp/ 10^) m'^, 
where t is the pulsar period time (in seconds) and Bp 
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is the surface pulsar magnetic field (in Tesla). The 
pair plasma density is expected to satisfy rig = rip = 
Mugj, where a moderate estimate of the multiplic- 
ity gives Af = 10 [2lj. Choosing this value and let- 
ting r = 1 s, we note that for the weak magnetar field 
strength Bp = 10^ T, the characteristic pump frequency 
w_Lchar ~ Wp/^, where the QED and plasma effects are 
of equal importance when the photon splitting process 
is of the order of o^^char ~ 4 x lO-'^^rad/s, i.e. in the 
optical range. For uj±_ <C LUA_chan the plasma nonlin- 
earities dominate, whereas the QED effect dominates in 
the opposite regime. The evolution of the coupled sys- 
tem of Eqs. (j4|) is. to a large extent, controlled by the 
pulse length of the pump mode. For long pulse-lengths 
with L ^ {uj±CE±/Eciit)^^ (Fig. 1), the system shows 
a "predator-prey" type of behavior where the energy os- 
cillates chaotically between the different modes. For a 
moderate pulse-length with L ^ {LJ±CE±/Ecrit)^^ , the 
excited EM wave energy propagates out of the interaction 
region, and one encounters a somewhat more ordered be- 
havior and an effective damping of the pump mode. The 
EM wave energy is then mainly converted to the paral- 
lel polarized forward scattered EM mode (Fig. 2). For 
a short pulse- length with L ^ {uj±CEj_/ Eait)"^ , ther- 
mal fluctuations do not grow due to convcctivc stabiliza- 
tion (the growth rate F vanishes for large wavenumbcrs), 
and effectively the nonlinear interaction vanishes. As the 
plasma density increases, plasma effects become impor- 
tant, as depicted in Fig. 3. For ujp ~ ^u;j_char, there are 
two simultaneous effects of the plasma. First, in this 
regime the QED and plasma contributions to the cou- 
pling strength are comparable, i.e. Cpi ^ Cqed, which in- 
creases the total coupling strength C = Cpi-l-CqED, since 
the phases of Cpi and Cqed coincide. Second, when the 
group velocity of the backscattered EM wave is slowed 
down, the effectiveness of convective stabilization is di- 
minished, increasing the interaction strength and speed- 
ing up the conversion of the EM wave energy. Simi- 
larly to the case without the plasma (Fig. 2), the EM 
wave energy mainly ends up in the parallel polarized for- 
ward scattered EM mode, but the characteristic splitting 
timescale is considerably faster with the plasma present. 
We stress that the simulations presented in this Letter 
contain results that are easily generalizable to other pa- 
rameter ranges, following the discussion presented above. 
The combined effect of plasma and QED effects should 
have consequences for the emission spectra from magne- 
tars and pulsars. While the QED effects alone shifts the 
spectrum towards linear polarization, we emphasize that 
the effect is much more pronounced when plasma effects 
are present, which holds for radiation with frequencies 
of the order tJchar ^ ^p/^j i-e. in the optical range for 
magnetars, while in the infrared to microwave range for 
pulsars. Thus, we suggest that evidence for a combined 
plasma-QED photon splitting process should be sought 
for in the polarization signature of magnetar and pulsar 



emission. 
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